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SUMMARY 


A method  for  measuring  in-plane  displacement  using  two  laser 
speckle  photographs  is  described.  The  technique  allows  the  user  to 
cancel  rigid  body  translation  and  rotation.  This  increases  the 
measurement  accuracy  of  relative  displacements  associated  with  load 
induced  strains. 

The  technique  involves  making  separate  speckle  photographs  of 
a test  specimen.  One  is  made  with  no  load  on  the  specimen  (no  dis- 
placement); the  second  is  made  with  the  specimen  loaded  (strained). 
A sandwich  is  constructed  from  the  two  speckle  photographs  and  data 
are  recovered  in  a manner  similar  to  that  used  with  conventional 
speckle  photography.  However,  the  fringes  generated  during  data 
recovery  are  not  strictly  Young’s  fringes.  Therefore  a curvature 
of  various  degrees  is  encountered.  A method  of  interpreting  the 
fringe  pattern  is  described  fully  in  the  text. 
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SECTION  I 


INTRODUCTION 


Laser  speckle  photography  has  emerged  in  recent  years  as  a profi- 
cient technique  for  measuring  small,  in-plane  displacements  (References 
1,  2,  and  3).  It  is  particularly  useful  for  determining  the  detailed 
characteristics  of  displacement  or  strain  fields  in  a small  area.  For 
example,  the  local  displacement  field  near  a tapered  shank  interference 
fit  fastener  has  been  measured  using  this  new  technique  (Reference  4). 

Accurate  determination  of  strain  from  displacement  field  measure- 
ments requires  better  than  average  precision  because  of  the  differen- 
tiations involved.  This  becomes  particularly  difficult  when  rigid  body 
displacements  are  present.  The  objective  of  this  report  is  to  present 
and  discuss  a new  method  of  using  speckle  photography  which  allows  for 
the  introduction  of  arbitrary  amounts  of  rigid  body  in-plane  displace- 
ment and  rotation.  This  provides  the  user  with  an  opportunity  to  cancel 
rigid  body  movements  of  a test  specimen  and  thereby  significantly 
increase  the  accuracy  of  relative  displacement  measurements  obtained 
within  a local  area. 

The  new  technique  is  called  ,?Dual  Plate  Speckle  Photography”.  The 
basic  idea  is  by  no  means  novel  and  has  probably  occurred  to  others. 
However,  there  are  some  fundamental  and  experimental  problems  involved 
which  will  be  discussed  in  the  text.  The  incentive  to  develop  this 
technique  was  provided  by  Dr.  Nils  Abramson  from  the  Royal  Institute 
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of  Technology,  Stockholm,  Sweden.  During  his  visit  to  the  Air  Force 
Flight  Dynamics  Laboratory  in  April  1974,  he  discussed  in  detail  the 
experimental  methods  he  developed  for  doing  "Sandwich  Holography" 
(Reference  5) . 

In  Section  II,  basic  elements  of  speckle  photography  are  briefly 
reviewed.  Two  of  the  appendixes  support  the  material  presented  in 
this  section  and  are  recommended  reading  for  those  not  previously 
exposed  to  these  ideas.  Dual  plate  speckle  photography  is  described 
in  Section  III.  A mathematical  model  is  developed  to  explain  the 
fringe  pattern  one  observes  in  the  diffraction  halo  of  a dual  plate 
specklegram.  Also  some  experimental  data  are  provided. 
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SECTION  II 


BASIC  SPECKLE  PHOTOGRAPHY 


The  speckled  appearance  of  a diffuse  surface  illuminated  by 
laser  light  is  easily  recorded  using  high  resolution  photographic 
film.  If  a double  exposure  is  made  of  an  object  with  some  in-plane 
displacement  occurring  between  exposures,  then  each  speckle  will  be 
recorded  twice.  The  distance  between  each  speckle  pair  is  a direct 
measurement  of  the  in-plane  vector  displacement.  The  photographic 
film  plate  on  which  such  a recording  is  made  is  usually  called  a 
"specklegram" . 

Displacement  data  are  recovered  from  a specklegram  by  passing 
a narrow  collimated  laser  beam  through  the  photographic  image.  High 
spatial  frequencies  in  the  speckle  pattern  cause  diffraction  and 
part  of  the  laser  beam  diverges  into  a cone  called  the  diffraction 
halo.  Since  the  speckles  occur  as  equally  spaced  pairs  in  an  other- 
wise random  pattern,  the  diffraction  halo  will  be  modulated  by  Young’s 
fringes.  These  fringes  will  have  a spacing  which  is  inversely  pro- 
portional to  the  speckle  pair  separation  and  perpendicular  to  it’s 
direction.  Thus,  displacement  magnitude  and  direction  can  be  measured 
at  any  location  on  the  image  by  passing  a narrow  collimated  laser  beam 
through  that  point.  Figure  1 illustrates  the  basic  recording  and  data 
readout  technique.  For  more  detailed  information  concerning  speckle 
photography,  the  reader  is  referred  to  Reference  3. 


The  generation  of  Young’s  fringes  in  the  diffraction  halo  is 
predicted  by  scalar  wave  theory.  The  calculation  is  closely  related 
to  a Fourier  transform.  A brief  sketch  of  a mathematical  model  is 
presented  in  Appendixes  A and  B. 
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SECTION  III 


DUAL  PLATE  SPECKLE  PHOTOGRAPHY 


In  order  to  cancel  the  effect  of  rigid  body  in-plane  movement  of 
a specimen  when  using  conventional  speckle  photography,  the  film  holder 
must  be  translated  or  rotated  between  exposures  by  a corresponding  amount. 
An  alternative  procedure  would  be  to  induce  an  opposite  but  equal  movement 
of  the  image  by  some  optical  means.  In  either  case,  prior  knowledge  of 
the  amount  of  rigid  body  movement  is  required. 

Dual  plate  speckle  photography  offers  the  user  a method  of  canceling 
rigid  body  movements  without  a prior  knowledge  of  the  magnitude.  Instead 
of  employing  a double  exposure  technique,  two  separate  speckle  photographs 
are  made.  The  film  plates  are  then  sandwiched  together  and  data  recovery 
is  accomplished  as  with  a conventional  specklegram.  In-plane  rigid  body 
translation  or  rotations  are  canceled  by  moving  one  film  plate  with 
respect  to  the  other.  This  allows  small  changes  in  relative  displacement 
as  a function  of  position  to  be  accurately  measured. 

With  a conventional  double  exposure  specklegram,  in-plane  displace- 
ments result  in  speckle  pairs  contained  in  a plane  perpendicular  to  the 
optical  axis.  This  geometry  is  not  completely  simulated  by  the  dual  plate 
technique.  Film  emulsions  are  not  infinitely  thin  and  glass  plates  vary 
in  thickness.  As  a result,  each  speckle  of  a pair  will  be  recorded  either 
forward  or  behind  itTs  partner.  This  is  true  even  when  the  recordings  are 
made  such  that  the  sandwiches  are  put  together  emulsion  to  emulsion.  This 
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situation  has  fundamental  implications  in  both  the  recording  and  recovery 
of  speckle  data.  Each  of  these  will  be  discussed  in  some  detail. 

The  following  procedure  is  employed  to  record  a dual  plate  speckle- 
gram.  A glass  film  plate  (we  typically  use  Agfa  10E75)  is  sandwiched 
with  a clear  plate  of  glass  placing  the  emulsion  side  towards  the  center 
of  the  sandwich.  A speckle  photograph  is  made  of  the  undisplaced  object 
with  the  film  plate  on  the  lens  side  of  the  sandwich.  Then,  second  speckle 
photograph  is  made  of  the  displaced  object  using  a second  identical  sand- 
wich except  that  it  is  turned  around  so  that  the  clear  glass  plate  is 
positioned  towards  the  lens.  Both  film  plates  are  then  processed  and 
thoroughly  dried.  A dual  plate  specklegram  is  constructed  by  placing  the 
two  speckle  photographs  emulsion  to  emulsion  such  that  the  images  are 
geometrically  matched.  A holding  device  to  accomplish  this  task  is  highly 
recommended.  Figure  2 is  a photograph  of  a dual  plate  holder  which  pro- 
vides for  translation  or  rotation  of  one  speckle  photograph  with  respect 
to  the  other.  A description  of  this  device  is  provided  in  Appendix  C. 

There  are  two  fundamental  limitations  in  the  recording  method  pre- 
viously described.  Both  are  related  to  the  fact  that  the  emulsions  on 
each  film  plate  may  be  in  slightly  different  planes  when  the  exposures 
are  made.  First,  laser  speckles  are  three-dimensional  entities  in  the 
image  plane.  The  two-dimensional  speckle  pattern  decorrelates  in  a dis- 
tance which  is  small  compared  to  the  conventional  depth  of  field  for  an 
imaging  lens  system  (References  2 and  6) . Decorrelation  is  particularly 
severe  with  lens  apertures  wider  than  about  f/4.  This  problem  can  be 
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reduced  by  increasing  image  magnification.  Second,  even  within  the 
speckle  correlation  distance,  a small  change  in  magnification  occurs 
with  image  position.  This  causes  a radial  shift  of  speckle  position 
for  image  points  off  of  the  optical  axis.  During  data  recovery,  this 
will  induce  apparent  in-plane  displacement  fringes.  This  phenomena 
is  closely  related  to  the  apparent  in-plane  displacement  fringes  that 
are  generated  with  a conventional  specklegram  when  out-of-plane  dis- 
placements of  the  object  have  occurred  between  exposures  (References 
2 and  6) . This  problem  can  be  minimized  by  using  imaging  lenses  with 
long  focal  lengths. 

Data  recovery  from  a dual  plate  specklegram  is  also  complicated  by 
the  longitudinal  displacement  of  paired  speckles  due  to  separate  emul- 
sions. The  fringes  modulating  the  diffraction  halo  are  found  to  be 
curved  in  various  degrees.  Typical  patterns  are  shown  in  Figure  3. 
These  photographs  suggest  that  a modification  of  the  Young’s  equations 
derived  in  Appendixes  A and  B are  required  in  order  to  interpret  these 
fringes  and  acquire  displacement  data. 

The  Fraunhofer  diffraction  integral  introduced  in  Appendix  A is 
restricted  to  apertures  contained  in  the  x-  <7  plane  at  Z=*o.  A physical 
interpretation  of  Eq.  (A-l)  is  that  it  represents  a summation  of 
Huygen’s  wave  amplitudes  at  an  arbitrary  observation  point.  The 
Huygen’s  wave  source  points  are  contained  in  the  aperture.  These 
are  temporally  coherent  since  it  is  assumed  that  the  aperture  is 
backlighted  with  a plane  wave  at  normal  incidence. 


Figure  3.  Typical  Fringe  Patterns 
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Now  interference  phenomena  in  the  diffraction  halo  are  primarily 
dependent  on  the  phase  summation  and,  therefore,  on  the  aperture  to 
observation  point  distance  given  by  Eq.  (A-2).  If  source  points  in 
the  aperture  occurs  for  z/o,  it  is  this  expression  that  requires 
modification.  Two  changes  are  necessary.  First,  the  source  point  to 
observation  point  distance  is  decreased  by  approximately  Z ( Cjo^l  O +* 
C&q.<$> -2.).  Second,  the  radiative  pha^e  of  the  source  points  are  now  a 
function  of  Z.  This  can  be  accounted  for  by  increasing  the  path  length 
by  a distance  equal  to  Z.  The  modified  expression  to  replace  Eq.  (A-2) 
is,  therefore, 

£ =■  50  — X ^<Usyt  & - y <f>  - z ( <p  - 2.  ) (1) 

One  further  change  is  required.  In  Appendixes  A and  B,  planar  apertures 
allow  summations  to  be  performed  with  an  area  element  c/xc/y  , With 
source  points  at  z^o,  the  area  element  can  be  a function  of  y and  Z. 
For  an  aperture  consisting  of  two  or  more  pinholes  represented  by  three- 
dimensional  Dirac  delta  functions,  a volume  element  dx  Jyclz  is  appro- 
priate. 

With  the  above  modification,  the  fringe  pattern  for  a dual  plate 
specklegram  is  easily  derived.  Let  a speckle  pair  be  represented  by  an 
aperture  W given  by 


^ = s(„ 


A- 


Z -h 


■tr)  + + £ , </,  z~  j;\ 


(2) 


11 


The  in-plane  displacement  A is  taken  in  the  X direction.  The  longitudinal 
separation  b models  the  situation  of  each  speckle  being  in  a different 
emulsion  plane. 

Following  the  procedure  outlined  in  Appendix  B - Case  I,  an  expres- 
sion for  the  intensity  J is  obtained. 

J oc  1 + Cscni  [ A k -sus*.  O -b  b k ( Z - ocn  G - o<^(t>)^\  (3) 

Examination  of  Eq.  (3)  reveals  a set  of  curved  fringes  of  the  type  shown 
in  Figure  3, 

Three  features  about  Eq.  (3)  should  be  noted  immediately.  First,  if 
b = 0,  Eq.  (3)  and  Eq.  (B-3)  are  identical  and  Youngfs  fringes  are  seen. 

If  A b this  situation  is  closely  approached  and  fringes  similar  to 
those  in  the  top  photograph  of  Figure  3 are  obtained.  Secondly,  if  A and 
b are  of  nearly  equal  magnitude,  the  pattern  is  substantially  altered. 

The  [bk  ocn.<t>3  term  causes  curvature;  the  term  modifies  the 

spacing.  A typical  pattern  is  shown  in  the  center  photograph  of  Figure  3. 
Finally,  if  A -C>  the  fringe  function  approximates  a set  of  concentric 
circles.  The  spacing  between  these  circles  is  inversely  proportional  to 
fc>.  The  bottom  photograph  in  Figure  3 illustrates  this  pattern. 

The  previous  discussion  indicates  that  although  data  recovery  is 
somewhat  complicated  by  the  longitudinal  separation  £>  of  paired  speckles, 
all  of  the  information  needed  to  measure  the  in-plane  displacement  A is 
present.  The  following  procedure  is  typical.  First,  the  value  of  b is 
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measured  by  translating  one  speckle  photograph  with  respect  to  the  other 
until  the  pattern  shown  in  the  lower  photograph  in  Figure  3 is  obtained. 
Fringe  minima  of  order  /i  are  then  related  to  b by  the  equation 

k = A (4 

f 1 - C<riL  © - Oe^L(p) 

Because  of  circular  symmetry,  Eq.  (4)  reduces  to 


(yi-  l/x}  A 

( 1 - Oc y-Q.  0 ) 


(5) 


Relative  values  of  in-plane  displacement  A can  now  be  acquired  at 
any  point  on  the  image.  The  value  of  A at  the  point  used  to  measure  b is 
taken  as  zero.  All  other  image  points  will  produce  fringe  patterns  with 
minima  of  order  )2  related  to  A by  the  equation 


A 


- b ( 3-~  ox-Q  - <f>  ) 


(6) 


If  measurements  are  secured  along  the  symmetry  line  defined  by  4>*^> 
Eq.  6 reduces  to 


A = 


(H-  ka-  ) \ — U 1 - CXWL  & ) 


(7) 


A study  of  Eqs.  (6)  and  (7)  reveals  several  subtle  properties.  Note, 
for  example,  that  for  positive  values  of  >7,  a negative  value  of  ©would 
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appear  to  yield  a negative  value  of  A.  This  is  not  a valid  interpreta- 
tion. Negative  values  of  9 correspond  to  negative  interference  orders.* 
The  sign  of  © can  be  determined  from  the  diffraction  halo  fringe  pattern. 
Positive  values  of  B sweep  the  observation  point  away  from  the  center  of 
curvature  (to  the  right  on  the  center  photograph  in  Figure  3) . 

The  sign  complications  discussed  previously  are  not  present  when  data 
are  read  from  a conventional  double  exposure  specklegram.  This  is  because 
Young’s  fringes  are  symmetrical,  and  therefore,  the  direction  of  displace- 
ment is  undeterminable.  Without  some  external  information,  it  is  not 
possible  to  know  which  speckle  in  a given  pair  was  recorded  first. 

With  dual  plate  speckle  photography,  displacement  direction  becomes 
determinable.  This  information  is  contained  in  the  curvature  of  the 
fringes.  If  the  front  plate  in  the  dual  plate  specklegram  is  a recording 
of  the  undisplaced  object,  the  relative  displacements  suffered  by  the 
image  recorded  on  the  rear  plate  are  in  a direction  pointing  away  from 
the  center  of  curvature  (positive  values  of  0). 


*Positvie  orders  are  denoted  by  71  = 1,  2,  3 •••  The  first  negative  order 
minimum  occurs  for  ft  = 0,  the  second  negative  minimum  occurs  for  71-  -1,  etc. 
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APPENDIX  A 


FRAUNHOFER  DIFFRACTION  INTEGRAL 


Consider  an  aperture  W of  finite  extent  illuminated  from  behind 
by  a collimated  laser  beam  (a  monochromatic  plane  parallel  light 
wave) . The  far  field  intensity  distribution  of  diffracted  light 
is  closely  related  to  the  Fourier  transform  integral  (Reference  7). 

In  particular  the  intensity  I is  proportional  to  a scalar  product  UlC 
where  U is  the  Fraunhofer  diffraction  integral. 


Equation  (A-l)  represents  a summation  over  an  aperture  W contained  in 
the  X-y  plane  at  2.-  o.  The  wave  number  k Is  £ /T  divided  by  the  wavelength 
and  S is  the  distance  from  a point  in  the  aperture  to  the  point  of  obser- 
vation. 

The  validity  of  Eq  (A-l)  is  restricted  to  small  angle  diffraction 
so  that  if  the  origin  is  located  near  the  center  of  the  aperture 


S ~ So  - x 0 - y l (f>  (A-2) 

where  $0  is  the  origin  to  observation  point  distance.  The  angles  O 
and  ({)  are  projections  into  the  x-z  and  y-*  planes  respectively  of 
the  inclination  of  the  ray  £o  from  the  x axis  (See  Figure  4).  The 
coordinates  X and  y locate  source  points  in  the  aperture  w.  Therefore, 
(A-l)  can  be  approximated  by 
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Figure  4.  Sketch  Defining  Angles  0 and  <j> 
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u 


-Ik  (*  ^ e 
JL 


(A- 3) 


Within  a multiplicative  constant  Eq  (A-3)  is  a two-dimensional  Fourier 
transform  of  the  aperture  W.  The  spatial  variables  are  x and  y.  The 
frequency  variables  are  i9  and  cux  - k st**.  <f>  , thus 


u 


= c 7f  WCx^)} 


(A- 4) 


where 


VJ 


-1  (o>,  x + 


(A-5) 


The  constant 


C = J-7T ' - 


JL 


JLkSe 
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APPENDIX  B 


YOUNG'S  FRINGES 


The  equations  developed  in  Appendix  A are  used  here  to  calculate 
the  intensity  distribution  for  two  different  apertures.  For  the  first 
case  the  aperture  is  two  small  pinholes.  The  resulting  diffraction 
pattern  is  shown  to  be  the  classical  Young's  fringes.  The  second 
calculation  models  the  aperture  as  a small  area  in  a double  exposure 
specklegram  illuminated  by  a collimated  laser  beam.  The  result  is  a 
random  diffraction  halo  modulated  by  Young's  fringes. 


Case  I - Two  Pinholes.  Consider  the  aperture  W to  be  two  very  small 
pinholes  separated  by  a distance  A.  This  aperture  can  be  modeled  by 
the  summation  of  two  Dirac  delta  functions  of  the  form  , y - ^ 

Thus 


IX  — C f x - *4, ) y ) ^ eT  ( x *-  A/x ) i ) y 

Evaluation  of  the  Fourier  transform  yields 

'LL  ~ X C Ocxl  (-JT  °*'  ) 

Now  the  intensity  J is  proportional  to  liVi thus 


or 


Jf  cC  1 -h  C*XL  ( A CO,) 

T °C  JZ  + O^-CL  ( A k & ) 


(B-l) 


(B-2) 


(B-3) 
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Equation  (B-3)  is  the  intensity  distribution  for  Young’s  Fringes.  By 
using  where  A is  the  wavelength  and  noting  that  minima  occur 

for  ( Ak  0)=— 1,  the  conditions  for  fringe  minima  are 

in.  - /a.  ) A = A AMVI  G (B-4) 

71  = l.ljl  • • • 

Case  II  - Specklegram.  Before  calculating  the  diffraction  pattern  from 
a double  exposure  specklegram,  first  consider  a single  exposure  photo- 
graphic transparency  of  an  object  illuminated  by  laser  light.  The 
recorded  image  is,  of  course,  speckled.  If  a small  diameter  laser  beam 
is  passed  through  a portion  of  the  transparency,  the  aperture  W will  be 
the  amplitude  transmission  function  inside  the  illuminated 

area.  Thus 

U0  = C 7[  Ru>i)}  (B-5) 

and 

R(x,  yj}  (B-6) 

Equations  (B-5)  and  (B-6)  can  be  evaluated  to  some  extent  since  the 
statistical  properties  of  the  speckle  pattern  y;  are  known 

(Reference  8).  This  will  not  be  done  because  the  objective  here  is 
simply  to  obtain  representative  expressions  for  the  diffraction  halo. 
Equations  (B-5)  and  (B-6)  are  these  expressions  and  are  denoted  by 
the  zero  subscript. 


X - ?{Ru,J  7* 
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Now,  for  a double  exposure  photographic  transparency  (speckle- 
gram)  where  the  object  has  suffered  an  in-plane  displacement  A in 
the  X direction  between  exposures,  the  amplitude  transmission  function 
is 


R (x-  7 ) ^ R(xi vA/'L)'l) 

A series  expansion  of  Eq  (B— 7 ) yields 


(B-7) 


where  primes  indicate  differentials  with  respect  to  the  variable  X. 
The  linear  properties  of  a Fourier  transform  provide  that 


*7 ($(*]}  -(*&>,) 


(B-9) 


(V) 

where  q ^ ^ is  the  7i  M,  derivative  of  a function  q (x)  • Thus 


VL 


} 


li 


C 3-l-Z 


H n 

» » • 

-j 


(B-10) 
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Using  Eqs  (B-5)  and  (B-6)  and  CO  t - k 


'LL  - UL  0 e2_  CsO-CL.  ( k.  /insist  & ) 


(B-ll) 


and 


I - Zc  j~  1 ■+■  Go-tf.  (a  s4As>s\  & ) j (B-12) 

Thus,  the  intensity  distribution  from  a specklegram  is  the  diffraction 
halo  from  a single  exposure  transparency  modulated  by  Young's  fringes. 
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APPENDIX  C 


DUAL  PLATE  HOLDER 


A dual  plate  specklegram  is  made  by  placing  two  speckle  photo- 
graphs emulsion  to  emulsion  such  that  the  images  are  geometrically 
matched*  In-plane  movement  of  one  film  plate  with  respect  to  the 
other  is  used  to  cancel  rigid  body  displacements  as  explained  in 
the  text.  A holding  device  to  accomplish  this  task  is  described 
in  this  appendix. 

A photograph  of  the  plate  holder  is  shown  in  Figure  2.  This 
device  has  three  main  parts;  a main  frame,  a translation  stage,  and 
a fixed  plate  holder.  Figure  5 is  a photograph  of  the  disassembled 
main  frame.  The  front  and  rear  sections  sandwich  a translation  stage 
and  a center  section.  Screw  adjustments  and  return  springs  are  housed 
in  the  center  section.  Figure  6 is  a photograph  of  the  translation 
stage  set  within  the  center  section.  A four  by  five  inch  glass  film 
plate  is  held  in  place  by  small  spring  clips  along  the  edge  of  the 
"window".  Figure  7 shows  the  fixed  plate  holder  removed  from  the  main 
frame.  It  too  is  equipped  with  spring  clips  to  hold  a glass  film  plate. 
This  fixed  plate  holder  was  designed  to  be  removed  so  as  to  gain  access 
to  the  translation  stage.  This  facilitates  mounting  and  removal  of  film 
plates  on  both  the  translation  stage  and  the  fixed  plate  holder.  The 
fixed  plate  holder  is  secured  to  the  main  frame  by  six  Teflon  cams. 

These  provide  adjustable  pressure  to  the  glass  film  plates  which  must 
slide  one  with  respect  to  the  other.  A light  silicon  lubricant  on  the 
film  emulsion  was  found  to  reduce  friction. 
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ge  Within  Center 


Figure  7 . Holding  Device  With  Fixed  Plate  Holder  Removed 
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